Efficient photon sorter in a high-dimensional Hilbert space 
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An increase in the dimension of Hilbert space for quantum key distribution (QKD) can decrease 
its fidelity requirements while also increasing its bandwidth. A significant obstacle for QKD with 
qudits (d > 3) has been an efficient and practical quantum state sorter for photons with complex 
wavefronts. We propose such a sorter based on a multiplexed thick hologram constructed from 
photo-thermal refractive glass. We validate this approach using coupled-mode theory to simulate a 
holographic sorter for states spanned by three planewaves. The utility of such a sorter for broader 
quantum information processing applications can be substantial. 
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We concern ourselves here with the secure distribution 
of a one time only key from a sender (Alice) to a re- 
ceiver (Bob). The three elements of any quantum key 
distribution (QKD) system are: (1) Alice must be able 
to prepare at will a single photon state chosen from a 
set of mutually- unbiased bases (MUB) [l|, 0], (2) each of 
these quantum amplitudes must be propagated from Al- 
ice to Bob with reasonable fidelity, and finally (3) Bob 
must have the ability to choose between one or another 
of the MUBs and, if he is lucky, be able to efficiently 
sort and detect each of these photon states. This QKD 
scenario has been exhaustively studied in the literature 
and is replete with security proofs for numerous proto- 
cols 0, U m ■ These security proofs have been extended 
in many cases to higher-dimensional Hilbert spaces Q, 
and all of the protocols have been or are currently being 
implemented successfully 0|. 

Conventional realizations of QKD today involve trans- 
mitting heavily-attenuated laser pulses from Alice to Bob 
and encoding qubit information in each packet by utiliz- 
ing the spin of the photon. This allows Alice and Bob, 
who are suitably authenticated, the possibility to estab- 
lish and share an arbitrarily-secure one-time only key be- 
tween them. Here they have access to a two-dimensional 
Hilbert space and can therefore form three MUBs each 
with two orthogonal polarization states. Such a six-state 
QKD scheme 3 has limited bandwidth and optical fi- 
delity constraints. These constraints can be ameliorated 
by extending the QKD to higher-dimensional Hilbert 
space 0. 

The potential of extending photon-based QKD to 
higher dimensions was made possible in 1992 when Allen 
et. al.0 showed that Laguerre-Gaussian light beams pos- 
sessed a quantized orbital angular momentum (QAM) of 
Ih per photon. This opened up an arbitrarily high di- 
mensional quantum space to a single photon[loj . Fol- 
lowing this discovery Mair et. al. [ll|, unequivocally 
demonstrated the quantum nature of photon 0AM by 
showing that pairs of 0AM photons can be entangled 
using parametric down conversion. Shortly thereafter. 



Molina- Terriza et. al. ^IS*! introduced a scheme to pre- 
pare photons in multidimensional vector states of 0AM 
commencing 0AM QKD. Recently a practical method 
has been demonstrated to produce arbitrary 0AM MUB 
states using computer-generated holography with a single 
spatial light modulator (SLM) [l3 |. 

While the advantage of 0AM QKD lies in its ability 
to increase bandwidth while simultaneously tolerating a 
higher bit error rate (BER) Q, two potential problems 
confront this approach. First, such transverse photon 
wave functions are more fragile in propagation than the 
photon's spin 15, [l^, and the divergence of the states 
i'xVl) may require larger apertures. Despite this, multi- 
conjugate adaptive optical communication channels may 
be able to ameliorate these problems [3 ■ A second 
obstacle involves the efficient sorting of 0AM MUB-state 
photons with small Fock-state quantum numbers, and 
the paper will address this particular problem. Currently, 
the only solution to this problem is the use if a cascaded 
Mach-Zehnder interferometric system 1^ 2^, 21 1. Pro- 
posed systems of this kind have been demonstrated only 
for 4-dimensions and are simply not practical. Other ap- 
proaches that use crossed thin diffraction gratings are not 
efficient enough to establish a secure key. 

We focus here on the efficient sorting of single pho- 
tons with arbitrary complex wavefronts. Ideally, what is 
needed is an 0AM version of a polarizer, i.e. a single op- 
tical element, one per MUB basis, that can efficiently sort 
each of the qudit states in that basis while equally dis- 
tributing every other qudit state. Thick holographic grat- 
ings fortunatel y p roduce high diffraction efficiency in the 

first order [13, Mm. If several predominant diffracted 

orders are required, as is the case for sorting, several inde- 
pendent fringe structures can exist in the emulsion. Such 
multiplexed holograms have been used for multiple-beam 
splitters and recombiners 24] and more recently for wide- 
angle beam steering |25j . In this paper we propose such 
a MUB-state sorter based on a multiplexed thick holo- 
graphic element constructed from commercially available 
photo thermal refractive (PTR) glass [2^. Due to the 
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unique properties of PTR glass the grating's thickness 
can approach several mm and be highly Bragg selective. 
There is evidence that such sorters can be highly effi- 
cient, > 95% [25|- Our simulations presented here and 



empirical data on thick Bragg gratings indicate that they 
may provide an adequate solution to this critical and long 
standing problem. 

Before we describe our proposed thick holographic 
MUB sorter we will briefly focus our attention on twelve- 
state QKD and work in a 3-dimensional Hilbert space. 
This is one more dimension than that available to pho- 
ton polarization states and serves to illustrate our ap- 
proach. Nevertheless, our work is equally applicable to 
higher dimensional Hilbert spaces. Its limitations will 
require further investigation. 

In 3-dimensions there are a maximum of four MUBs 
which we refer to here as MUBi, MUB2, MUB3 and 
MUB4. Each of these orthonormal bases contain three 
state vectors. If we identify \a), \b) and |c) as the or- 
thonormal ket vectors of MUBi, then the other nine 
qudit states from the other three MUB bases are specific 
linear combinations of these (Table |T| . 

For our application we can freely choose as our AIU Bi 
any three pure 0AM states (|a), \b) and |c)) correspond- 
ing to an angular momentum, la — ah, h = bh, and 
Ic — ch with integers a, b and c being the azimuthal 
quantum numbers. For the purpose of our calculations 
we can simplify our analysis and retain the physical con- 
tent by quantizing in the space of linear momentum (k- 
QKD) rather than in angular momentum (OAM-QKD). 
As a result we can freely choose as MUBi any three 
non co-linear planewaves. In this case, our three integer 
quantum numbers will be the number of waves of tilt of 
these planewaves with respect to the normal of the holo- 
graphic emulsion of aperture D. These waves correspond 
to a transverse linear momentum = ahk^ , = bhk^ 
and pjf = chk^\ respectively. Here, = k\/D is the 
x-component of a plane wave of frequency with one wave 
of tilt (r ^ A/I?). In the frame of the hologram and in 
units where the speed of light is unity, the components 
of the 4- momentum {p — {p* ,p^ ,py ,p^}) of each of our 
three photons can be expressed in terms of their trans- 
verse momentum k^ and wavenumber (k). 



Pa ^hka ^ ahk-'il, 1,0, ^{k/ak^y 

Pb =hkb = bhk^{l, 1, 0, ^/{k/bk^f - 1}, 

^hkc =cftA:^{l,l,0, \/(fcM^)^-l}- 

In the remainder of this paper the transverse linear mo- 
mentum wavenumbers represent our three quantum num- 
bers for k-QKD. These three planewaves define our first 
MUB. 

MUBi^{\a),\b),\c)}. 

Each of these states represents a transverse Fourier mode 
of a photon; they are orthogonal — Sij) and define 




FIG. 1: An illustration of our proposed thick fiolographic 
MUBi sorter. The three signal waves are the appropriate 
linear superpositions of the three qudit states of MU B4, shown 
in Table m 



our 3-dimensional Hilbert space. The other nine MUB 
states (Table H]) can be obtained from these by linear su- 
perposition and will represent wavefronts with both am- 
plitude and phase variations. 

For each MUB we consider a multiplexed thick holo- 
graphic sorter, i. e. a triple-exposed grating structure 
formed by the incoherent superposition of three gratings 
within a single emulsion. Here, each grating is formed 
by the superposition of the respective MUB state and its 
own unique plane reference wave. In this paper we con- 
centrate on the construction of the MU B4 sorter (Fig.[T]), 
as the other three MUB sorters will be of similar design. 

We construct the MUB4 sorting in three steps. First, 
we record the interference pattern of our first signal wave 
[04) with a corresponding reference planewave, |ri) hav- 
ing ri ^ 1 waves of tilt. After this initial recording is 
complete, we then record the interference pattern of our 
second signal state from MUB 4, namely I64) with a sec- 
ond reference planewave (|r2). Finally, we record a third 
independent set of fringe patterns by interfering the qudit 
signal state |c4) with a third reference planewave jra). 
This produces a triple-multiplexed hologram. 

We show that the hologram described above faithfully 
represents a quantum projection operator for MUB4. 



7^4 = |ri)(a4| + k2)(64| + k3)(c4| 



(1) 



Its operation on any one of the 12 MUB states should 
produce the desired result. In other words, if the holo- 
gram is illuminated by MUB state [04), I64) or |c4) it 
should produce a planewave in state \ri), \r2) or jra), re- 
spectively. If it is illuminated by any of the other nine 
MUB states it should then produce an equally weighted 



3 



MUBi 


MUB2 


MUBi 


MUBi 


la) 
\b) 
\c) 


\a2) oc |a)+ 16)+ \c) 
62) oc \a) +z 16) \c) 
\c2) OC la) + z^ 6) + z c) 


as) oc a)+ 6) + 2 c) 
63) oc a) + 2 6)+ c) 
IC3) oc a) + 2^ 6) + 2^ \c) 


04) oc la)+ 16) + 2^ |c) 
64) oc la) + 2 6) + 2 |c) 
C4) oc |a) + 2^ 6)+ Ic) 



TABLE I: The four MUBs in our 3-dimensional Hilbert space. Note that phase 2 = exp(i27r/3) is a cube root of unity, and we 
have suppressed the normalizing factor of l/v^S in each of the nine MUB states in the last three columns. 



response into all three reference states, e.g. 

{r,\Vi\h)^^, V*G {1,2,3}. 

The unique property of PTR glass with its bulk in- 
dex of refraction, no — 1.4865, and depth of modulation, 
An/riQ = 336ppTO, place it squarely in the realms of 
scalar diffraction theory and coupled-mode (CM) theory. 
Furthermore, a PTR hologram can be thick D cm 
with Bragg-plane periods, A ~ A. Consequently, their 
Bragg selectivity ^ A/L can approach the diffraction 
limit of one wave of tilt across its aperture [2^. For 
A = 1085 nm and Z? ~ L ~ 1 cm this yields a minimal 
divergence of our three signal waves of a few arcsec. 

To examine the MU B4 sorter we will follow closely the 
CM approach of Kogelnik ^23l| and the notation used by 
Case [24] to solve the scalar wave equation for polariza- 
tion perpendicular to the plane of incidence. 







(2) 



Here, the linearly-polarized electric field Ey{x,z) of fre- 
quency v is assumed to be independent of y. Following 
CM theory we keep only primary modes for the electric 
field. These are the transverse harmonic modes given 
by the k- vectors ki, fc2, fcs, ka, ki, and kc associated to 
planewave reference states |r2), jra) and MUBi sig- 
nal states I a), |c), respectively. 



Ey{x,z) =i?i(z)exp' 



''^■^+R2{z) e^p''^-^ +R3iz) exp'=3-r 
+Saiz) exp''--^ +Sbiz) exp''"-^ +Sciz) exp^-■^^ 



Here, the six mode amplitudes, {Ri} and {Si} are only 
functions of z. They are set initially to {Ri{z = 0)} — 
{1, 1, 1} and to the corresponding amplitude and phase 
factors of one of the twelve corresponding signal states 
shown in Tabic [D For example, signal state |c4 > 
one would set {Si} — l/-\/3{l, z^, 1}. The wavenumber 
k{x, z) of Eq.[2]represents the three incoherently recorded 
gratings mentioned above. 



n{x, z)kQ 



noko^ 

13 



An {Iri + Ir2 + Im) 
no 6(1 + 



where Im is the intensity modulation of the i*"^ grating, 
e.g. 



/fl3 = 2 - 



+ ze 



ikh-r 



Assuming that the interaction between the diffracted or- 
ders is slow, we can neglect second order terms and arrive 
at the CM equations for the mode amplitudes. 
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where 



= (- 



-), and 



kf and cr, = fcf 



V6(3+V3)/ 

are the z-components of the wave vectors for the refer- 
ence and signal states, respectively. The solution of this 
equation for each of the twelve initial signal states (MUB 
states) are shown in Fig. O and faithfully reproduce the 
desired projection operator of Eq. [TJ We independently 
examined the far-field pattern for such gratings using a 
finite difference time domain solution of Maxwell's equa- 
tions and observed that the CW assumptions were valid, 
i.e. only the primary modes were dominant, and there 
were no relevant polarization changes in the field. 

While the analysis presented here suggests that a high 
efficiency single optical element sorter is feasible with 
commercially available materials and holographic record- 
ing techniques, further work is needed. First, we need 
to understand the quantum nature of such an element. 
Second, we need to examine its sensitivity to alignment 
(linear and rotational) and examine the wavelength scal- 
ing issues involved in recording 0AM photons. Our first 
step will be to produce a 3-state k-QKD MUBi sorter. 
Our goal is to test its performance using states generated 
by a single phase SLM 1^. If the MUB-state sorters 
described here can be produced, they should have far 
more utility in quantum information processing than just 
QKD, e.g. as an essential element in linear quantum com- 
puting [27]. 
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FIG. 2: Here we show the predicted performance using CM theory for the multiplexed MUB4 sorter. This single optical element 
is a triple multiplexed volume hologram constructed using PTR glass (no — 1.4865, An = 0.0005 at A = 1085 nm). This small 
ratio of depth of modulation to bulk index {An /no ~ 336ppm) places it squarely in the linear regime of CM theory. There is one 
graph for each of the twelve MUB states which show the probability that the signal photon will be observed to be diffracted into 
one or another of the three original reference directions {{r\ri) — Rie^''^'^ , {'r\r2) — i?2e**^ ' and (rl^s) = R^e^^^"^) as a function 
of the depth {z)oi the emulsion. The curves have been terminated near their maximum efficiency depth {zmax ~ 8.5 mm). 
Each of the three curves in each graph has its own unique value of Zmax- They can be made roughly equivalent by proper 
choice of wavevectors p^'s and the cr^'s. The solid curve corresponds to \r\), the dashed curve to |r2) and the dotted curve to 
Ira). A large probability amplitude in the |ri) corresponds to a signal wave with a large projection along the qudit state \ai). 
The same is true for the other two MUB4 states. The twelve graphs are portrayed in a matrix. Each column represents each 
of the four allowable MUB bases, and each row corresponds to one of the three MUB states within each MUB basis (Table |T]|. 
For this MUB4, sorter we see perfect sorting in the last column yielding the desired quantum state projection. However for any 
photon prepared in any one of the qudit states in the first three columns, there is equal probability that the photon will be 
diffracted into the three directions (]ri), |r2) and jrs)) thereby yielding no information as to the identity of the quantum state. 
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